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We study the role that both vacuum fluctuations and vacuum entanglement of a scalar held 
play in identifying the spacetime topology, which is not prescribed from hrst principles—neither in 
general relativity or quantum gravity. We analyze how the entanglement and observable correlations 
acquired between two particle detectors are sensitive to the spatial topology of spacetime. We 
examine the detector’s time evolution to all orders in perturbation theory and then study the 
phenomenon of vacuum entanglement harvesting in Minkowski spacetime and two hat topologically 
distinct spacetimes constructed from identihcations of the Minkowski space. We show that, for 
instance, if the spatial topology induces a preferred direction, this direction may be inferred from 
the dependence of correlations between the two detectors on their orientation. We therefore show 
that vacuum huctuations and vacuum entanglement harvesting makes it, in principle, possible to 
distinguish spacetimes with identical local geometry that differ only in their topology. 

PACS numbers: 04.62.+V, 04.20.Gz, 03.65.Ud 


I. INTRODUCTION 

The field equations of general relativity are local. As 
such, they describe the local structure of spacetimes, 
while remaining silent about the large scale structure 
of our Universe, including its spatial topology [T]. Dif¬ 
ferent cosmological models produce a variety of global 
properties, including a range of spacetime topologies. In 
principle, we do not expect a future theory of quantum 
gravity to fix the spatial topology of the Universe. For 
example, loop quantum gravity admits the possibility of 
topological changes [5] , while path integrals in spin-foam 
models are summed over all possible topologies [5]- How 
our actual observable Universe behaves, from the point 
of view of topology may be derivable from the initial con¬ 
ditions of the Universe [3] , but at present we believe that 
the underlying fundamental theories are compatible with 
whatever large scale structure we happen to live in. 

It is then a task of observational astrophysics and cos¬ 
mology to determine the actually realized scenario 
Unfortunately, to date, astronomical observations of the 
large scale structure of the Universe are compatible with 
either an Euclidian spatial topology or any kind of open 
or closed topology where the relevant topological scales 
are comparable or above the Hubble radius. Even if the 
signals we are collecting may have originated in causally 
disconnected regions, our entire measurement record is 
overwhelmingly local in space and time. 

Quantum mechanics possesses inherently non-local 
properties. The paramount manifestation of those prop¬ 
erties is embodied by quantum entanglement [5117] . The 
fundamental question that we address in this paper is: 
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Can we use quantum mechanics to amplify the non¬ 
locality of our experiments and thus become more sensi¬ 
tive to aspects such as the topology of spacetime? 

From the point of view of local observers the vacuum 
state of any quantum field is entangled, and thus localized 
vacuum fluctuations are correlated [5]. While the Unruh 
effect EHin] is perhaps the most spectacular manifesta¬ 
tion of vacuum entanglement, acceleration of observers 
is not required. It was demonstrated that correlations 
in the vacuum measured by local inertial observers are, 
in principle, strong enough to violate Bell-type inequal¬ 
ities m- Furthermore, it is indeed known that local¬ 
ized particle detectors can extract entanglement form the 
vacuum state of a quantum field, even while remaining 
spacelike separated, due to the intrinsic non-locality of 
the field ground state [IMS]. This phenomenon has 
become known as ‘entanglement harvesting ’ m, and it 
has already been proven to encode information about the 
global structure of spacetime. Specifically, Ver Steeg and 
Menicucci demonstrated that entanglement harvest¬ 
ing is sensitive to the large scale structure of spacetime. 

Continuing the study of the dependence of vacuum en¬ 
tanglement of a quantum field on cosmological parame¬ 
ters, we consider a real massless scalar field and a pair 
of two-level Unruh-DeWitt detectors [m uni and de¬ 
rive their density matrix assuming only existence of the 
Wightman function. Using Minkowski spacetime as a 
benchmark, we compare the resulting entanglement be¬ 
tween the two detectors, with the counterparts in topo¬ 
logically non-trivial locally flat spacetimes. 

This paper is organized as follows: In Section [T^ we 
review the detector model and derive the joint density 
matrix of two detectors in Minkowski space, initially in 
their ground state and at rest with respect to one an¬ 
other, to all orders of perturbation theory and present 
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a consistent leading-order expression. We then analyze 
the probability that a single detector undergoes a transi¬ 
tion to the excited state, and quantify the entanglement 
harvested between the two detectors resulting from inter¬ 
acting with the field. In Section [lll| we consider an anal¬ 
ogous situation in two locally flat, topologically distinct, 
quotient spacetimes constructed from Minkowski space. 
Since entanglement measures are locally unobservable, 
we evaluate and compare, in addition to local statistics, 
the observable correlations between the detectors. We 
conclude in Section m with a discussion of the results 
presented and outline future directions of research. 

We use the H-convention for the signature of the 

metric and work in natural units h = c = 1. 


II. DETECTORS IN MINKOWSKI SPACE 

In this section we analyze general properties of two 
Unruh-DeWitt detectors A and B that interact locally 
with a real massless scalar field for (effectively) a fi¬ 
nite period of time. We assume that the spacetime in 
question allows us to separate solutions of the Klein- 
Gordon equation into positive- and negative-frequency 
modes, thus allowing us to define the Wightman function 
[9], and that the two detectors are situated in the same 
coordinate patch. This is enough to obtain the struc¬ 
ture of their joint detector density matrix p^B (Eq- (|^)j 
and a general form of its matrix elements and entangle¬ 
ment properties. Specific information about investigated 
spacetimes enables us to calculate this expressions ex¬ 
plicitly. While our analysis is based on perturbation the¬ 
ory, we obtain results that are valid at all orders. As a 
benchmark we provide the leading order expressions for 
two identical detectors at rest in Minkowski space. 


A. The detector model 

We will model our particle detectors as two-level quan¬ 
tum systems that interact with the real massless scalar 
field (p via an Unruh-DeWitt monopole coupling [18]. Un¬ 
der this coupling model, the time-dependent interaction 
Hamiltonian in the interaction picture is given by 

Hi (t) = A(r) (e*^^cr+ -f e"*^^cr") (p [a;(r)], (1) 

where r is the proper time of the detector, A (r) is a 
weak time-dependent coupling parameter which controls 
the strength and length of the interaction, U is the energy 
gap between the ground state |0)^ of the detector and its 


excited state |1)^, (t=*= are SU(2) ladder operators which 
act on the state of the detector, that is <7+ |0)^ = 11)^;, 
a~ |1)^ = |0)^, (cr=*=)^ = 0 and (P{x{t)) is the field eval¬ 
uated along the trajectory of the detector. This model, 
while simple, captures most of the relevant features of 
the light-matter interaction when no angular momentum 
exchange is involved [20l |2T] . 

As we will assume the detectors are at rest with respect 
to one another, it is convenient to parameterize both de¬ 
tectors and the field’s time evolution by the common co¬ 
ordinate time t. We express the coupling parameter as 
X(t) = eocpt), where eg 0 is the coupling strength and 
e{t) = is a Gaussian switching function. Due to 

the strong suppression for |t| ^ a, we can approximate 
the detector being “on” when \t\ < a and “off” otherwise. 

Prior to the interaction, we consider that the detec¬ 
tors are initially in their ground states |0)^ and |0)^, 
and the field is in the vacuum state |0), so that the ini¬ 
tially joint state of the two detectors and field is given by 
1^) “ |0 )a |0)b |0)- During the interaction the composite 
system undergoes the unitary evolution 

JJ — J'g-i I 

where T denotes time ordering and the Hamiltonians Ha 
and Hb describe the field interaction with detectors A 
and B, respectively. These Hamiltonians are given by 
Eq. 0- In principle, the detectors may have different in¬ 
teraction coupling strengths, switching functions, and en¬ 
ergy gaps. Even when considering identical detectors we 
will use the subscripts A and B for the coupling parame¬ 
ters = cb = eg, since it allows us to easily distinguish 
local and non-local terms in the perturbation expansion 
below. The two interaction Hamiltonians commute at 
equal times, [HA{t), HBit)] = 0. 

Some of the calculations are more conveniently per¬ 
formed using the Schrodinger picture. We label the basis 
state of the detectors as mA,mB S {0,1} and the states 
of the field as p. The combined system is given by 

I'*/’(^)) = CmArriB 

mA,mB,p 

e-<^^A+E-'-B)*\mAmB)U^{t)\p) , (3) 

where the (generalized) basis states |^) of the scalar field 
are freely evolved by a unitary operator and the sum¬ 
mation and integration with appropriate measures are 
performed over the entire multimode Fock space. Hence 
the index p can be 0 (the vacuum state), or take a single- 
or multi-mode labels. The coefficients Ck^kBi^it) the 
solutions of 
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iC-k^kB vif) — 


^ ^ {I^ I </> [xA{t)] I /i) 


+ ^Cfc^,fcB+iMWe^ (!/|0[a;B(i)] I A^) 


(4) 


where we use mod 2 addition for the detector indices kA 
and ks- 

We will assume the two detectors are initially unex¬ 
cited and the field to be in the vacuum state, which cor¬ 
responds to the initial state |^') = |0)^ |0)^ |0) . This 
initial state corresponds to a single non-zero coefficient 
Cqqq = 1 at the zeroth order in the perturbative expan¬ 
sion of the coefficients. Beginning with the zeroth order, 
higher order corrections to the coefficients CkAkBv(t) are 
obtained by substituting coefficients of one lower order on 
the right hand side of Eq. In what follows we take 
t —> oo, which yields the coefficients in the late time limit 
when the interaction between the detectors and held has 
ceased. 

The joint state of the two detectors is pab = 
tr^jU ('k| W], where the trace is over the held degrees 
of freedom. It follows from Eq Q that the nth order coef- 
hcient where p involves m < n particles, depends 

on the coefficients of the order n — I with the opposite 
detector parity sum kA + ks and mil particles. As 
a result, at all orders of perturbation theory the density 
matrix 

{PAB)kAlB,rnAnB ' ^^aIb aub u (®) 

has the form of the so-called X-state [25], 


Indeed, tracing out either of the detectors in the state 
PAB, say detector B, results in the state pA of detector 
A 


PA 


1-A 0 
0 A 


(9) 


in the basis {|0)^,|1)^}. The same result is obtained 
when we analyze a single detector that interacts with 
the held via the Hamiltonian IIa- If the appropriately 
dehned separation between the detectors (e.g. distance L 
between the stationary detectors) increases, they become 
virtually independent and the total state approaches the 
direct product of the density matrices of the individual 
detectors, 


/I-A 0\ ^ n-B 0\ 
P.4B-(^0 b) ■ 

As a result, AT —)• 0, C —> 0, and B -A AB. 

To simplify the exposition we consider the case of two 
identical detectors, i.e., Ba = Hb- Hence, in translation- 
invariant spacetimes we substitute H > A in Eq. ([^. In 
the leading order of the perturbation parameters tA,^B, 
the quantities A, X, C and E are given by 


Pab = 


/ r-ii 

0 

0 

ri4e“*«\ 


A 

0 

r22 

^■236“*^ 

0 

(6) 


0 

?’23e*^ 

r22 

0 


\ri4e*« 

0 

0 

r44, / 




' —OO J —<50 




in the basis {|00), |01), |10) , |11)} where \ij) = \i) A\j) b, 
and all the coefficients are positive. Since pab 
is a valid density matrix, the normalization condition 
= 1, and the following two positivity conditions 
must be satished: 


X = —2eAtB / dt / dt' e{t)e{t') 

J —OO J —<50 

B^i^+^’)W{XA {t),XB {t')) + 0{et), 


C = CA^B 


dt / dt' 


( 11 ) 


( 12 ) 


riiTA^ > rf4, r22r33 > r^^. 


(7) 


A useful form of this matrix that explicitly separates 
the local and nonlocal quantities is 


PAB = 



AO 0 A:^ 

B-E C 0 

C* A-E 0 

0 0 Ej 


( 8 ) 


XBit)) +0{et) , (13) 

E=\X\^ + A^ + 2\C\^+0{€^o), (14) 

where W {x,x') := (0| 0 (x) (/) (x') |0) is the Wightman 
function. For an explicit derivation of pab to all orders in 
perturbation theory for the general case when EIa ^ Hb 
see Appendix fA| 

In four-dimensional Minkowski space the Wightman 
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function for a real massless scalar field is given by 


Wm (x, x') = -—: sgn S({x — x')^) 

47rj 

1 


{x — x'Y 


(15) 


where {x — x'Y = VtJ.vix''' — x''^){x'' — x'"). Following 
[231 - 125] we treat the Wightman function as a distribu¬ 
tion, regularizing the resulting integrals via the principal 
value prescription as necessary; see Appendix]^ for more 
details. 

We consider the two detectors to be at rest with respect 
to one another and separated by a distance L = |x^—| ■ 
We find the matrix elements of pab to be 


^2 r 


Am = 


e erfc (an) -|-C>(eo), (16) 


47r - 


Am = 

dyTT L 


l+erf|jF) 


Oief 


0 / > 
(17) 


Cm = 


4y/w L 


— e 4 <t 


Im 


gifiL erf \ i 2 —^ 


-sm{nL) +0{e^), (18) 


has a negative eigenvalue. This is known as the Peres- 
Horodecki criterion |7| . A popular measure of the entan¬ 
glement between the two detectors, commonly used in 
the context of entanglement harvesting, is the negativity 

11—1 

JC (pab) ■= - - 2 - •='-P-^ ( 21 ) 

where || • || denotes the trace norm, and p- is the negative 
eigenvalue in question [7]. 

Another measure of entanglement that we employ is 
the concurrence, which plays an important role in entan¬ 
glement theory and allows us to quantify the amount of 
entanglement in terms of the maximally entangled pairs 
of qubits (ebits) [7|. For an arbitrary two-qubit state p 
the concurrence is equal to 

C(p) = max( 0 , Ai - A 2 - As - A 4 ), ( 22 ) 


where the A^’s are the square roots of the eigenvalues of 
the matrix pp, where p = {ay ® ay)p*{ay (g) ay) with ay 
being the Pauli y matrix, ordered such that Ai > A 2 > 
A 3 > A 4 [55]. We discuss the relationship between these 


two measures of entanglement in Appendix A 2 


From the Peres-Horodecki criterion it follows that X- 
states are entangled if and only if either of the alterna¬ 
tives 


Em — I Am I + -I- 2(7^ -I- C>{eY), (19) 


tIa > ^ 22 X 23 , ?’23 > ^ll'’44, (23) 


where erf( 2 :) is the error function, erfc(z) = 1 — erf(z), 
and the subscript M denotes that these quantities are 
calculated in Minkowski space. We outline the details 
of this calculation in Appendix A1 and B. In the limit 
of infinite interaction time, i.e., when the two inertial 
detectors are always on, the resulting transition rate is 
zero as expected [26]. In the limit of large L the density 
matrix indeed approaches the direct product pa® Pb- 


holds. 

For two identical detectors in the state paB given in 
Eq. ([^, these conditions are equivalent to 


\X\-A + O{et)>0, \C\-^ + O{4)>0, (24) 


respectively. However, Eqs. (18) and (nQ)) ensure that the 


second condition is never satisfied. Hence the entangle¬ 
ment is 


B. Information-theoretical properties of the joint 
state 

Prior to the interaction with the field, the two detectors 
and the field were initially in a separable state jdi) = 
|0)^ |0)^ |0). After the interaction the detectors become 
entangled. Since each detector interacts locally with the 
field, and there is no direct detector-detector interaction, 
any resulting entanglement between the two detectors 
must have been redistributed from entanglement already 
present in the vacuum state of the field. This is the 
phenomena of entanglement harvesting [l 2 l [T3l [T5l ITHl 
m], and it can be thought of as a kind of entanglement 
‘swapping’ from the vacuum state of the field to the joint 
state of the detectors. 

A two-qubit state p is entangled if and only if its par¬ 
tially transposed matrix p ^^, 


c/2 =A7 = max(0,|A| - A-hC>(e;))) , (25) 


if and only if ri 4 > r 22 - Using the expressions in Eqs. 
(161 and ([T^ , we plot the concurrence of the state pab 


in Fig. |la| 

For a system of two qubits, the concurrence can be used 
to calculate the entanglement of formation Ep, defined 
as the number of the maximally entangled states needed 
to prepare pab [7] 


Ep {pab) = h 


1 + y^l — C { pab )^ 




41n2 


(l-ln(CV4)) +0(. 


oj ^ 


(26) 


where h{x) = —a;log 2 a; — (1 — 1 ) log 2 (l — x). From the 
above we see that in the standard units of a bipartite en¬ 
tanglement, the so called ebit, the entanglement between 
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the two detectors in a perturbative regime is relatively 
week: it scales as Eq. 

Neither of these quantities is accessible by local mea¬ 
surements of either detector. Instead, we consider two 
local observers, one in possession of detector A and the 
other in possession of detector B, each measuring the ob¬ 
servable Pauli z operator , and quantify the correlation 
in the outcomes of their measurements. We characterize 
their results by random variables ta and r^, respectively, 
withr^, vb G {Oj !}• The correlation between these vari¬ 
ables is given by 


cov^B 

corr^B = - 

ctactb 


E- AB 

"^^A{l-A)B{l-By 


(27) 


where coy ab '■= {i'ai’b) — (?'a) (?'b) is the covariance be¬ 
tween TA and tb and <j\ = coy a a and (7^ = covbb are 
the variances associated with ta and tb- In Minkowski 
space, when the switching function of the two detectors 
is coincident, the correlation between the measurements 
performed on the two detectors is given by 

corr" =li5£4±22lL + 0(.3). ( 28 ) 


III. DETECTORS IN QUOTIENTS OF 
MINKOWSKI SPACE 

We now consider an analogous situation in two differ¬ 
ent quotient spacetimes of Minkowski space with non¬ 
trivial spatial topologies. The first is a cylindrical uni¬ 
verse Mo 


Mo = M/To, (29) 

which is constructed as a quotient of Minkowski space 
M with the group Tq = {J//} generated by the discrete 
isometry Jq : (t, x, y, z) i—>■ (t, x,y, z + £). 

The second spacetime we consider M-, is again a 
cylindrical universe in which rotations by tt in the xy- 
plane have been identified 

M-=M/r_, (30) 

where the group r_ = {J"} is generated by the discrete 
isometry J_ : {t, x, y, z) i—>■ {t, —x, —y, z -\- £). 

The compactification scale of these two spacetimes is £. 
Both Jo and J_ preserve space and time orientation and 
act freely and properly Q this ensures that both Mo and 
M- are space and time orientable Lorentzian manifolds. 


^ Let G be a group which acts on a topological space X. A group 
action G X X X is said to act freely if for all x & X, gx = x 
implies g is the identity. A group action G X X X X X is 
said to act properly if the map {g, x) i—)■ {gx, x) is proper, that is 
inverses of compact sets are compact. 



0 12 3 4 

(a) The concurrence C (pab) /^q detectors in Minkowski 
space 


3.0 


2.5 



0 12 3 4 

aQ 

(b) The correlation function corr^g/cQ for detectors in 
Minkowski space 

FIG. 1. The (a) concurrence C (pab) /cq and the (b) corre¬ 
lation function corr/is/co are plotted for two identical detec¬ 
tors in Minkowski space A4, as a function of their separation 
L and their energy gap Q in units of a. The red line in (b) 
corresponds to = Am, above which the concurrence 

vanishes exactly. The solid black regions in the plots indi¬ 
cate the plot has been clipped at (a) corr ab/eq = 2 and (b) 
Cm {pab) /cq ~ 0.25. 


As neither Jq and J_ affect the Minkowski line element, 
both Mo and M- are locally flat spacetimes. The be¬ 
haviour of particle detectors in both of these spacetimes 
has been studied in the past, with a focus on applications 
to Hawking radiation and the Unruh effect [iniES]- 
The advantage of considering spacetime topologies 
built from quotients of Minkowski space, is that the 
Wightman function in the quotient spacetimes can be 
easily constructed via the method of images from the 
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Wightman function in Minkowski space. If we are 
given a Green’s function Gm{x,x') on a spacetime A^, 
we can construct the corresponding Green’s function 
Gm/j{x,x') on the quotient spacetime AijJ by the im¬ 
age sum m 


oo 


Gm/j{X:x')= ^ rrGM{x,rx'), 


n——oo 


(31) 


where Jx' denotes the group action of the group element 
J on x' and 77 = {1,-1} corresponding to normal scalar 
fields and twisted fields respectively. From now on we 
will restrict ourselves to 77 = 1 for simplicity; the results 
presented are easily generalizable to twisted fields. 

Making use of Eq. along with the Wightman func¬ 
tion in Minkowski space, Eq. (15), the Wightman func¬ 


tion in both J\4o and M- can be constructed, and used to 
evaluate the transition probability A in Eq. ( |11[), and the 
matrix elements X and C in Eqs. (12) and(|13|) respec¬ 


tively, in both spacetimes. We consider two detectors 
with the world lines 


XA{t) = {t,dA,ZA), XB{t) = {t,dB,ZB), (32) 

where and are two dimensional vectors lying in 
the xy-plane. Matrix elements of pab ap to the second 
order in cq in AIq are given by 




M ' 


X Im 


4i/7r 


n^O 


ni 


-inm f ^ ^ 

2cr 




— sin(nl!n) ) (33) 


An = A 


M 


4^/n 


Go — Cm + 


Im 


76 






e 


n^O 


1 -I- erf [ i — 

2(j 


(34) 


di/TT 


E 


Lr, 




erf 


-e 


2(7 


ail 


— sin {ilLn) 


Oiel) 


(35) 


where + rrr — 2n£Az and Az = za — zb\ while 


in A4_ these quantities are given by Eqs. (33)-(35), with 
the substitutions 


:=f"+4-|P(n), (36) 

L„ ^ Ll := + ri^£^ - 2nlAz + dd^^ 5 ^( 77 ), (37) 


where dk = ||dfc|| with k G {A,B} and P{n) is zero or 
one for even and odd n respectively. 

Note that A_, A_, and C- for detectors in A4_ de¬ 
pend on the absolute position of the two detectors, as can 


be seen from the dependence of Eqs. (36) and (37) on 


dyi and d^. This is a qualitative difference between A4_, 
and both the cylindrical universe A4o and Minkowski 
space Al, stemming from the fact that A4_ is not trans- 
lationally invariant in the xy-plane, which can bee seen 
from the isometry J_ used to define Al_. Gonsequently, 
A- ^ B-, as defined in Eq. ([^, which was true for both 
A1 and A4o, and Eq. (27) must be used to calculate the 
correlation function. 

From the Eqs. (33)-(35l we see for a large compactifi- 
cation scale £ of the quotient spacetimes, i.e., £ —>■ 00 , 
the contribution to the transition probability and ex¬ 
change probability from the image sums vanish, resulting 
in Aq = A- = Am and Xq = A_ = Xm- Gonsequently, 
in this limit the two universes are indistinguishable by 
measurements of either the transition probability of a sin¬ 
gle detector or the resulting entanglement between two 
detectors. 

To examine the effect the spatial topology of a space- 
time has on the transition probability of a detector, in 
Fig.i we plot the transition probabilities Am, Aq, and 
A_, in Minkowski space A4, the cylindrical spacetime 
A4o, and the spacetime A4_. The oscillations in the 
cases of A4o and A4_ are expected and akin to the ap¬ 
pearance of modified quasi-normal modes in spacetimes 
with closed topologies (see, for instance [5^1. 

The dependence of the transition rate, that is the 
derivative of the transition probability, on spatial topol¬ 
ogy for detectors that have been on for an infinite amount 
of time in both A4o and A4_ has been studied in [3D]. In 
black hole spacetimes with identical geometry but differ¬ 
ing spatial topology, the transition rate of detectors has 
been studied in [33] . 

In examining the effect the spatial topology has on 
vacuum entanglement, we plot in Fig. [^ the difference in 
correlations functions between detectors in A4 and A4o 
and detectors in A4 and AI _. 

In both Alo and A4_ the concurrence of the joint state 
of the two detectors depends on the orientation of the 
detectors with respect to the identified direction. Thus, 
in principle, by measuring vacuum corrections one may 
infer a preferred direction induced by the spatial topology 
of the universe. To illustrate this point, in Fig. [^ we 
plot the dependence of the concurrence of pab on the 
orientation of the two detectors in the AIq spacetime 
with respect to the identified direction. 

For the purposes of plotting, in Figs. m we truncate 
the image sums at n = ±10, as the inclusion of more 
terms does not affect the plots. This is because the con¬ 
tribution from larger n terms in the image sums decease 
quickly as n grows: oc e“" /n, which can be seen from 
Eqs. ([33|)-(|^. 


IV. DISCUSSION AND OUTLOOK 

In the spacetimes we analyzed (A4o and A4_), the ef¬ 
fects of global structure show up as small deviations in 
transition rate of a single particle detector, and entan- 
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— Minkowski — f^ja — 0.3 — ijo — 0.5 — ija — 1 


(a) Transition probability in yV)o 



— Minkowski — t/cr — 0.3 — tja — 0.5 — ^/cr = 1 


(b) Transition probability in Jv[_ for d^ja = 1 

FIG. 2. The transition probability of a detector in Minkowski 
space M is compared to the transition probability of a detec¬ 
tor in both (a) the AIo spacetime and (b) the M- spacetime, 
by plotting it as a function of the energy gap of the detector, 
in units of ct, for different circumferences I of the universe. In 
both (a) and (b) the dashed blue line is the Minkowski tran¬ 
sition probability Am- In the leading order, the probability 
of the transition |1) —>■ |0), i?i — i ?2 = 12 > 0, is equivalent to 
the excitation probability |0) —>■ |1) with < 0. The inset is 
a magnihcation of the excitation probability fl > 0. 


glement and observable correlations between two detec¬ 
tors from their counterparts in the spatially Euclidean 
Minkowski spacetime. However, in the limit of zero ex¬ 
trinsic curvature, i.e. when the compactification scale 
becomes infinite, these deviations approach zero and the 
results coincide with that of Minkowski space. We plan 
to disentangle the role of extrinsic curvature from that of 
entirely topological effects in future work. 

The Minkowskian result is recovered in the limit Ija ^ 
00 . Since in cosmological scenarios we expect the topo¬ 
logical scale to be at the order of the Hubble scale. 



0 12 3 4 


(a) Difference in corvAB between detectors in Ai and Ado 



0 12 3 4 

crfl 

(b) Difference in corr^^ between detectors in A4 and A4- 


FIG. 3. As a measure of the effect the spatial topology has 
on the observable correlations in pab, we plot the difference 
between the correlation function for detectors in Minkowski 
space and (a) the A4o spacetime: corr^^/eg — corr'^g/cQ, and 
(b) the A4- spacetime: corr^^a/cg — corr^g/eg- In addition, 
we plot the contour lines of the Minkowski correlation func¬ 
tion corr^^g/eg to serve as a benchmark. In both plots we 
have chosen both detectors to lie in the aiy-plane, so Az = 0, 
and choose the length of the identified direction to be ^ = 1. 
In plot (b) we have chosen dx and ds to be parallel so that 
L = dg — d-A, and have set dA = 0.1. The red region indi¬ 
cates when the correlations between the detectors is greater in 
Minkowski space than the identihed spacetime {A4o or M-) 
and vice versa for the blue region. The solid blue and red re¬ 
gions indicate where the plots have been clipped respectively 
at 0.10 and —0.10. 


































FIG. 4. The concurrence of pab for two detectors in the 
cylindrical nniverse Mo with circumference £ = 1 is plot¬ 
ted as a function of their orientation 6 with respect to the 
identified ^-direction. As Mo is translational invariant, we 
can choose our coordinate system such that detector A is 
at the spatial origin xa = 0 and detector B is located at 
xs = (L cos 0,1/sin 0). 

£ ^ ^/H and the effective experimental run-time is at 
the order of years at best, we expect tja > 10^°, making 
the effects unobservable. By the same token we expect 
that entanglement between different degrees of freedom 
during early Universe, where the relevant dimensionless 
ratios were of order of 1 or higher, to be significantly 
impacted by the emerging global structures. We plan to 
investigate these effects in future work. 

On the other hand, although the discussion throughout 
this paper has focused on the spatial topology of the en¬ 
tire Universe, the tools used and results presented apply 
equally well to fields and detectors in cavities with ap¬ 
propriate boundary conditions. Specifically, the results 
given regarding the cylindrical universe Mq are equiva¬ 
lent to detectors and the field in a cavity with periodic 
boundary conditions. Further, the blue regions in Fig. 
suggest that entanglement harvesting from a quantum 
field may be increased by constructing cavities with ap¬ 
propriate boundary conditions. In fact, there is already 
substantial evidence that this is the case: Entanglement 
harvesting in cavities has already been analyzed non- 
perturbatively in |34j and further in |35j . where it was 
shown that a combination of harvesting in cavity setups 
complemented with communication yields a sustainable 
source of quantum entanglement. This particular ampli¬ 
fication of harvesting in cavity setups is otherwise impos¬ 
sible in free space. 

A different group of questions deals with entanglement 
and correlations between the detectors in relative motion 
[5], as well as with the effects of delay between switch¬ 
ing the detectors. Finally, it is interesting to investigate 
the build-up of the correlations between the detectors in 
time. 


The density matrix of Eq. Q has the same form at 
all orders of perturbation theory. Nevertheless, it will be 
instructive to obtain non-perturbative results that are 
based on non-perturbative methods following the several 
different formalisms developed for harmonic oscillator de¬ 
tectors [aiiiMlETi, as well is to compare different types 
of detectors. 
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Appendix A: State of the detectors and 
entanglement measures 

1. Calculation of pAB 

Here we provide the details of the perturbative calcu¬ 
lations that lead to the reduced density matrix in Eq. 
Q. We present the derivation in Minkowski space, but 
it translates verbatim to any stationary spacetime. 

Here we work in the Scrbdinger picture as it will be 
more convenient than the interaction picture to extract 
the structure of the reduced density matrix pab to all or¬ 
ders in perturbation theory. The Hamiltonian describing 
two detectors interacting with the scalar field (j) is given 
by 

H = Ha + Hb+H^ + (Al) 

where Ha — and Hb = ^^bO'z are the free 

Hamiltonians of detectors A and B with energy gaps 
Ha and is the free Hamiltonian of the scalar 

field, and iJint is the interaction Hamiltonian in the 
Schordinger picture 

= eA{t)ax4>s{xA{t)) + eB{t)a^(j)s{xB{t)). (A2) 

where the subscript S reminds us that we are working 
in the Scrodinger picture. We choose the switching func¬ 
tions eA{t) and eB(t) to be proportional to our perturba¬ 
tion parameter eg 1; in the analysis presented we took 
eA{t) = e-Bit) = , but in general other switching 

functions may be studied. Since both detectors are in¬ 
ertial and at rest with respect to one another, we have 
equated the proper time of each detector with the coor¬ 
dinate time t. 

Initially (tin —oo) the detectors are in their ground 
state and the field in the vacuum state 

I'l') = |0 )a |0)b |0). (A3) 
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We expand the state IV'CO) = exp{—iHt) |'I') of the combined system in the eigenstates of the unperturbed Hamil¬ 
tonian at time t as 


\i^{t))= ^ (crnAniBtJ.{t)e 

rriAmB ^ 


(A4) 


where rrii S {0,1} labels the energy levels of the detectors with energies Ex^ — and Ei^ — Eq^ = ^b, 

the index /i represents a decomposition over the basis states of the scalar field, and U{t) = 0 Is 0 U^{t) where 

U^{t) = eyi'p{—iH^t). The coefficients CmAmB a{t) satisfy the Schrodinger equation 


'^mArriB H 


^ ^’'^'>*eAit)iaa:)mAnA (M I </> [XAit)] \ v) Cuatub ^ (0 ) 

TLA,I' 

+ i (/X | (/) (t)] k) ^ (t) ) , 


(A5) 


where ax is the Pauli x matrix. We represent the coefficients CmAmBU ^ more explicit form as 


= ( 


m.A'fyiB M — \^mn v i ^mnpi ^mn pip 2 i * ■ v 5 


(A6) 


i.e., spell out explicitly the vacuum, one-particle, two-particle, etc. components. To simplify the notation we label 
the terms Cgip as b{p) and cup^p^ as x{pi,p 2 ), and suppress the arguments. Contributions of these elements to the 
inner product that is performed with the appropriate measure will be denoted in the usual vector form, such as |a;p 
or b ■ a*. 


The solution to Eq. (A5), subject to the initial condition given in Eq. Eq. (A3), is given schematically as 


cooa — (l + -I-0, + C’(eo); Oj ^>(^ 0 )) (A7) 

com = ( 0 ,eo 6 W +e36(3),0,O(e3)) +0(e5), (A8) 

cioa = (Ojeoa*-^^ -I-0) C>(co)) +^(co)) (^9) 

cim = 0, + ©(eg). (AlO) 


As a result, the non-zero matrix elements of the reduced density matrix pab in Eq. ([^ describing the joint state 
of the two detectors, up to the fourth order in eg are given by 


Poo, 00 =: rii = 1 -f cosy + eg + 2gl‘^^ cos 74 -f | 

Poi,oi =: r 22 = eo|a^^^|^ + eo(a^^^ ’ 


Pro ,10 =: r 33 = -f • b^^^* + b^^^* ■ 5^^^), 


Pll,ll —: ^44 — Cg 


Poi,io =: r 23 e -I- eg(6*'^^ • -I- • b^^^*), 

Poo ,11 =: ri 4 e“*^ = 


-ZC = ,2.,(2),-*X2 + g4 (g(2)42)g*(7.-X2) + ^2) . ^(2)*^ _ 


(All) 

(A12) 

(A13) 

(A14) 

(A15) 

(A16) 


By construction pab is normalized. 

To calculate the transition probability, correlations, and measures of entanglement in the leading order, we need 
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the coefficients c 


■mAmBfj. up to order e§: 

/ OO 

dte^^^^esit) {p\(l){xB{t)) | O) , 

-OO 

/ OO 

{p\(l){xAit)) |0), 

-OO 


-- 


dt / dt' 


' —OO o —OO 


„( 2 ) 


Cl 1 ,, — 


/ dt dt' 

' —OO j —OO 


eAit)eAit')e~^^^^*~*'^ {v I (l){xAit))(l){xA{t')) I 0) 

+ ^ | 4>{xB{t))(l){xB(t')) \ O) 

+ | 4>{xBit))(l){xA{t')) \ O) 


(A17) 

(A18) 


(A19) 


(A20) 


where in the second order coefficients the index b stands either for vacuum (0) or a two-particle state (piP 2 )- Pio.iOj 
Poi.oij and poo.ii 

The matrix elements of pab are obtained from Eqs. (A17)-(A20) with the help of Eq. ([^. For identical detectors 
the matrix elements pio.io = Poi.oi) Poi.iOj and poo,ii are obtained straight forwardly from Eqs. (A17)-(A20), and 
result in the quantities A, X, and C respectively, and are given in Eqs. (11|-(13). 

As for the matrix element applying Eq. ([^ results in 


Pii.ii — dt dt' 


ft" 


dt" 




dt"'€{t)€{t')e{t")e{t"')e 


' —OO j —OO 


(O \(j){xB{t'"))4’{xA{t"))(l){xA{t))(j){xB{t')) \ O) -I- (O \(j){xA{t"'))4>{xB{t"))(t){xA{t))(j){xB{t')) \ O) 

+ (O \4>{xB{t"'))(j){xA(t"))(j){xB{t))(j){xA{t')) I O) + (O \(j){xA{t"'))4>{xB{t"))(l>{xB{t))(j){xA{t')) \ O) 


■ (A21) 


The four 4-point functions appearing in Eq. (A21) can be expanded in terms of products of Wightman functions 
using the commutation properties of the field, yielding 


nOO pt 

Pii.ii = / dt j dt' j dt 

J —cx> J —OO 


‘OO pt" 

ff I 


/ —OO J —OO 


-b 


dt'''e{t)e{t')t{t'')e{t''')e'^^^+^'-^''-^'''^I^W{xA{t'''),x 

W(xA(t''),XA(t))W(xB(t"'),XB(t')) + W(xA(t"'), XA(t))W(xB(t"), XB(t')) 

+ W(xA(t''),XA(t'))W(xB(t'"),XB(t)) +W(xA(t"'),XA(t'))W(xB(t''),XB(t)) 

(xB(t"'),XA(t))W(xA(t"),XB(t')) + W(xB(t"),XA(t))W(xA(t'"), XB(t')) 


(A22) 


where we have exploited the fact that the detectors are 
at rest with respect to one another, which results in 
W(xA(t),XB(t')) = W(xB(t),XA(t')). 


The first term appearing in Eq. (A22) is immediately 


identified as |A| , where X is given in Eq. (12). By 
changing the integration variables to 


,/// J. 

u = t — t, 

V, = t -ft, 


v = t" - t', 
v = t'' + t', 


(A23) 


and then to 


a = V — u 


b = V + u, 


(A24) 


the second and third terms in Eq. (A221 result in 
AB and 2C^ respectively, where A and B are defined 
in Eq. ( [TI| ) and C in Eq. ( |I^ . Thus we find 


pii.ii = \X\^+ AB + 2C‘^. 


(A25) 
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The matrix element poo^oo is obtained either directly 
from its definition or by exploiting the normalization con¬ 
dition iri^pAB) = 1 - 


2. Quantifying the entanglement in pAB 

We give here negativity and concurrence for the density 
matrix of (|^. In the case > r 22 r 23 

A/" = - 5 (r22 + ^33 - \/(r22 - r33)2 -I- 4ri42), (A26) 

and 

C = 2(ri4 — •y/r22C23), (A27) 


Eq. (B2| is not manifestly covariant, and requires addi¬ 
tional manipulations to obtain, for example, the Unruh 
effect [3S]. On the other hand, a straightforward cal¬ 
culation demonstrates that the distributional form of the 
Wightman function as given by Eq. (151 is free from these 
complications. 

For convenience, we summarize here several properties 
of distributions that we employed in obtaining Eqs. (16) 
and 0 m- Recall that the definition of a distribution 
G acting on a test function / is given by 


{G,f) 


9{y)fiy) dy, 


(B3) 


which simplifies further if r 22 = ^23 to 
Af = C/2 = ri 4 - r 22 - 
If r |3 > riir 44 then 

Af = -{rii + r44 - +r44Y +4r23), 


(A28) 


where the function g(y) defines the distribution G. The 
derivative of a distribution is obtained from the above 
definition by integrating by parts to give 


(A29) 


(G',/)=-(G,/'). (B4) 


and 


C = 2 (r23 - . (A30) 

In a general two-qubit system concurrence and nega¬ 
tivity are related by [35] 

C ^ 2AA > V(1 - Cy +C^ + {l-C). (A31) 

The negativity is equal to the concurrence if the eigenvec¬ 
tor of the partially transposed state corresponding 
to its negative eigenvalue is one of the Bell states (up 
to local unitary transformations). Indeed, for the identi¬ 
cal detectors (r 22 = ^ 23 ) when ri 4 > r 22 the two quan¬ 
tities coincide and the eigenvector in question becomes 
^(0,-1,1,0)^. 


Appendix B: The Wightman function 


To ease comparison between different sources, we first 
spell out how our convention for the Wightman function 
W := {0\(l){x)(j){y)\0) relates to the definitions in other 
references we use. In particular, 

W{x,x') = G'^{x,x') = —iD~(x,x') (Bl) 


where G'*' is introduced in [33] (and is called in the 
massless case), D~ is introduced in [33], and is called 

in [35] . 

Wightman functions are well-defined distributions, i.e. 
they can be represented as a distributional limit of reg¬ 
ular analytic functions [23] . However, representation in 
terms of functions that are also covariant requires a reg¬ 
ularization procedure, e.g., the Pauli-Villars regulariza¬ 
tion. A simple popular representation uses the fe” 
prescription [^ 


We{x — x') 


1 1 

dTT^ {t — t' — ieY ~ ~ 


(B2) 


The distribution 1/x acting on a test function fix) is 
defined as 


1 

5 

X 


fix) 



fix) 

5 

X 


(B5) 


where PV denotes that the principle value of the integral 
should be taken. All the subsequent inverse power dis¬ 
tributions l/x" are defined as distributional derivatives 
of 1/x, hence 

(/2./w) = (h/'w) 

^ r (B6) 

Eq. ( |B 6 [) is used in arriving at the expression for Xm 
given in (|17|). 

Particular care is required for the evaluation of inte¬ 
grals involving the delta-function term of the Wightman 
function when y := x — x' = 0. We consider the distri¬ 
bution 5(y^) as a limit of 

S (y^) = lim 6 (y^ - r^) 

= lim ;^((5(y-r)-f (5(y-hr)'). (B7) 

1—>0 2 r V / 


This prescription results in the action of the distribution 
sgn(y)i5(y^) on a test function /(y) to be 


/ OO 

dy sgn (y) 6 {y^) f (y) = f (0). (B 8 ) 

-00 


Eq. (B 8 ) is used in arriving at the expression for 
given in ( 161 . 
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